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Abstract 



We consider static spherically symmetric solutions of a general scalar 
field theory with non-standard kinetic energy coupled to gravity with a view 
to explain dark matter halos as a coherent state of the scalar field. Con- 
sistent solutions are found with a smooth scalar profile which can describe 
r"! \ observed rotation curves. A characteristic feature is that the rotation curves 

are linear for small r. The nontrivial solutions have negative scalar energy 
^SJ . density near the origin though the total energy is positive definite. We 

^ \ also reconsider the no scalar hair theorems for black holes with emphasis 

on asymptotic boundary conditions and super-luminal propagation. Some 
■rj" ■ modifications and extensions of previous analyses are discussed. 
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1 Introduction 

Recently there has been a lot of interest in the applications of non-canonical scalar field 
theories in cosmology. This is due mainly to the fact that the energy momentum tensor in 
such field theories has the potential to describe cosmological fluids with negative pressure 
which is a necessary ingredient for accelerated expansion. Examples are k-essence [T] which 
attempts to explain the accelerated expansion of the universe as well as the coincidence 
problem, the DBI action [2], tachyonic matter [3],|1], the ghost condensate model [5] and the 
Chaplygin gas model [6]. They are also interesting in the context of inflationary models [T,'8]. 
Since these theories contain non-standard kinetic energy terms, the constraints of global 
hyperbolicity and absence of super-luminal propagation play an important role P, (TUl [H] . 

In this paper we consider static, spherically symmetric solutions of the Einstein equations 
coupled to a general k-essence scalar field. The physical motivation is to look for consistent 
solutions describing galactic halos |T2] and black holes. The standard assumption concerning 
galactic dark matter halos is that they consist of an incoherent collection of weakly interacting 
massive particles. There have been attempts [131 [iH] where some very special type of scalar 
theories were used to discuss the possibility that the galactic halos could in fact be considered 
a coherent excitation of a scalar field, much like a boson star. This scenario would have the 
advantage of providing a unified treatment of dark matter and dark energy since the latter 
can be described by a k-essence like theory to begin with. 

We consider the most general scalar field lagrangian which depends on the field and 
invariants of the kind X = g^^d^(f)dy(f) formed from its first derivatives. We do not assume 
that the kinetic terms are separable nor that they are a quadratic form in the first derivatives. 
We find that solutions do exist which can describe galactic halos and for certain choices of 
the metric function give a good description of the observed rotation curves. Two classes of 
solutions are found: those that have negative energy density near the origin and those that 
don't. These can be phenomenologically distinguished by the shape of the rotation curves 
near the origin. Strictly speaking, only one of these solutions has a valid classical description, 
whereas the other may depend on new physics at short distance scales. For the case when the 
scalar energy density becomes negative near the origin, the total energy or mass enclosed can 
still be positive. Recent evidence from string theory indicates that there is no justification 
for restricting to potentials that are positive everywhere as long as the total energy is finite 
|14] . Negative energy densities are not unphysical by themselves as the Scharnhorst effect 
[15] clearly shows. Here one has faster than light propagation in a Casimir vacuum. In 
the literature there are studies [16j associating negative energy densities with super-luminal 
propagation of signals. However, in [11] it was argued that super-luminal propagation does 
not necessarily imply the violation of causality for which one requires closed time-like curves. 
All this clearly deserves further study to see if super-luminal propagation without causality 
violation is consistent. If so, then k-essence like theories would provide a larger context in 
which to study the formation of galactic halos. 

Another question discussed in this paper is the existence of black hole solutions in the 
combined gravity, k-essence system. This problem is of physical interest since galaxies are 
thought to have black holes at their center. Hence, we would like to ascertain whether black 



holes can coexist with the halos. Two possibihties are usually mentioned in connection with 
the issue of black hole-halo coexistence [TB] • The first is that the massive black holes were 
formed together with the galaxies through some internal dynamical process, or secondly that 
the black holes are primordial. In the second possibility, they were present even before any 
luminous activity, and in fact are the source driving the quasars. In both of these situations, 
once formed, the black holes continue to grow in time. If there is a scalar field pervading the 
universe then it can interact with the black holes and if the scalar no hair theorems [19J for 
static spherically symmetric black holes are valid it could cause their accretion. One could 
study this question by considering the stability of the halo in the presence of a background 
black hole solution. In this paper, however, we do not do this. Instead our more modest 
approach is to revisit the scalar no hair theorems. These are sometimes used to argue [18] 
that the black holes become heavier in time by "eating" the scalar hair. We present new 
ways to understand the theorems and suggest some avenues how they may be circumvented 
in the context of k-essence like theories. 

The paper is organized as follows. In section 2 we discuss the various steps that are 
necessary to model a spherically symmetric scalar halo and to describe the rotation curves. 
In section 3 we address the question of the no hair theorems for black holes. We conclude 
in section 4 with a discussion of our results. 



2 Scalar Fields and Dark Matter Halos 

In this section we will describe a halo assuming that it is made up only of scalar fields. 
Including an exponential disc of baryonic matter should not change the essential conclusions. 
Thus, we are interested in a scalar field theory coupled to gravity and described by an action 
of the general form. 



^ = y rf^Xv^L(X,0), (1) 

X = g^''V^(j)V.(j). (2) 

Its energy momentum tensor is given by, 

9 r C 

T^, = ^=^— = 2LxV^4>V,<i) - g,,,L. (3) 
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Here and elsewhere in this paper, Lx denotes the partial derivative with respect to X. The 
equation of motion for the scalar field is given by, 

G"^"V^V,0=^(L^-2XLx<^) (4) 

&'' = Lxg"" + 2LxxV^0V>. (5) 

The quantities satisfy several constraints which we now outline: an examination of the 
characteristics of the scalar equation of motion gives the speed of propagation of the scalar 
fluctuations, and demanding that this so called sound speed (C^) is not super-luminal imposes 



the constraint ^^^ > 0. Demanding in addition that the initial value problem be well posed, 
and the scalar field equation of motion be globally hyperbolic gives the following list of 
constraints for this system pfTUt [TT]. 

Lx > 0, Lxx > 0, Lx + 2XLxx > 0. (6) 

In addition, stability requires [TT] that Cs > 0. 



As discussed in the introduction we will be interested in static spherically symmetric so- 
lutions of the combined Einstein-scalar system of equations. These solutions must match 
on to the cosmological solution at large r. Depending on the model under consideration, 
the cosmological solution could either be almost asymptotically fiat or not. We will mostly 
assume almost asymptotically fiat solutions, but in the next section we will also consider 
other possibilities as well. The space-time metric will be described by, 

ds^ = e^dt" - e^dr^ - r" {dO^ + sin^edcj)^), (7) 

where z/ and A are functions of r alone. The (tt), (rr) and (66) components of the Einstein 
equations G(^ = nTj^ are given by: 

^{-l + e^ + rX') = -KL, (8) 






e^-ru') = -K{L-2XLx), (9) 

^ (rv + \\v'f + z/' - A' - irz/'A') = -kL. (10) 



2r ' 2 ^ ' 2 

The equation of motion for the scalar field (j4]) in the static limit may be written as, 

^" + [- + (- + -)( — )W + 1{ Vz^^:^} = 0, (11) 

^ ^^2^^2^ r^^Lx^2XLxx^ 2^Lx^2XLxx^ ^ ' 

where ip' = e~^(f)', and the speed of sound Cg is given by. 

Note that if we require the absence of super-luminal propagation of the scalar field fiuctua- 
tions then for all r, 6{r) > 1, and stability necessitates Cg > as well. 

We are now ready to discuss the question of galactic halos. In ref [20j, it is shown how some 
essential features of the metric function u can be deduced directly from the observed galactic 
rotation curves independent of the matter content and a specific gravitation lagrangian. From 
stability considerations [2D] concludes that. 



< ru'/2 < 1. (13) 

Further assuming circular halos and that information travels to us along null geodesies, one 
get's an additional constraint [20], which for our purposes may best be written as: 

\<(r)\ < ^^^. (14) 

^ = vl (15) 



In the above, Vc is essentially the tangential component of the velocity [21], i.e., the rotation 
curve. These inequalities indicate a smooth scalar profile. It is clear from these that it is 
not at all difficult to find functions Vc{r) which give a realistic representation of the rotation 
curves: one that grows for small r and flattens out in an intermediate region before matching 
on to the cosmological solution outside. This is precisely the kind of behavior we demand of 
the corresponding derivative of the metric function v so that our model is able to describe 
the observed rotation curves. 

Let us first look for consistent solutions of the different equations of motion near the origin. 
First, consider the scalar field equation (TTTj) . Since e^^ 7^ near the origin, we rewrite it 
in terms of the variable rj' = e~2 0'. The advantage is that now the scalar equation may be 
written without explicit reference to the metric function A: 

V" + -(^ + 2)<5(r)r7' + hir) = 0, (16) 

r 2 

We would like to determine the appropriate boundary conditions for small r. Since we are 
describing a halo as a coherent state of a scalar field, much like a boson star, the appropriate 
boundary condition at the origin is the generalized "no force" condition, i.e. ^ = 0, where 
p is the pressure and ^ is given by the Oppenheimer-Volkov equation. For our model, the 
pressure and the density are respectively, p = —L and p = L — 2XLx- Using (TTT]) or (TT6|) . 
it is straightforward to obtain: 

dp 2 ,rv' , ,9 , ^ 

-f = —{— + 2W'Lx. 18 

ar r 2 

Since the quantity inside the brackets is bounded between 2 and 3, we obtain the boundary 
condition for the scalar field equation that r]' goes to zero faster than ^/r at the origin (keep 
in mind that Lx > everywhere). We would next like to determine the exact rate at which 
this happens, so for this purpose consider a small r expansion for rj and write, 

r] = r'{co + Cir + C2r'^^ ). (19) 

As a consequence of the condition that rj' vanish at the origin, it follows that h{r) must 
have a smooth behavior at small r. To see why this is the case, we will first assume the 
opposite, i.e. suppose that h{r) goes like ^ to leading order at small r. The equation of 
motion flT6l) implies that to leading order rj' approaches a nonzero constant at small r, which 
rules out this possibility. Similarly, if /i ~ l/r", where n > 2 is an integer, then r]' ~ l/r""^, 
which would contradict our previous finding that 77' — ;> at the origin. So the only allowed 
behavior of h{r) consistent with the boundary conditions is that it is well defined at the 
origin. Substituting (TT9|) into (TT6l) . we find for the indicial equation (cq 7^ 0): 

s{s-l + ^5) = (20) 

2 < 7 < 3. (21) 

In the above, 7 denotes the value of 7 = rz/'/2 + 2 for small r and the inequality above 
follows from (IT3ll . Further, if we demand the absence of super-luminal propagation of the 



scalar fluctuations then S, the value of S for small r, is greater than 1. Thus we see that the 
absence of super-luminal propagation of the scalar fluctuations forces upon us the solution 
with s = 0. The actual behavior of t]' with r near the origin depends on the small r behavior 
of /i(r)|j. By substituting the expansion flT9l) with s = into flT6|) we find ci = if h(r) is well 
behaved near the origin, and it is non-zero when h{r) ~ 1/r. Thus depending on the small 
r behavior of h{r), Vj either vanishes at the origin or approaches a constant. In either case 
since Cq 7^ 0, r^ approaches a constant at small r. The boundary condition for our case is that 
77' — > at small r, thus the only behavior of /i(r) consistent with the absence of super-luminal 
propagation of the scalar fluctuations is that it is either a constant or vanishing at small r. 
In either case, C\ = 0, which is what we henceforth take to be the case. We will see below 
how these conditions translate into the shape of the rotation curves for small r. 

We must now find a consistent solution for small r for the metric functions as well. For 
these we turn to the Einstein equations (JHl) and ([9]). It is useful to rewrite these as 

A' = -KrL- (e^-l)(KrL + -), (22) 

z/' = nrp + (e^ - l)(/trp + -). (23) 

On the right hand side of these equations we have the energy density p = —L and the 
pressure p = L — 2XLx- From these we see that the combination (z/' + A') only depends on 
p — L, which is proportional to X (see also (13T1) below), u and A themselves depend on p and 
L separately. In general, the small r behaviors of L and {p — L) ~ X can be different: L is a 
function of both X and 0, and (j) goes to a constant at small r. From a Taylor expansion of L 
around X = Cg, it is seen that the leading behavior of L can therefore be either a constant or 
that of at least X. From fl22l) and fl23l) we see that in the former case, the leading behavior 
of z/' and A' are determined from that of L and this leading behavior is cancelled from the 
sum. It is the subleading terms of u' and A' which are now proportional to X. For small r, 
we have two possibilities: (i) both L and 2XLx have similar behavior near the origin, i.e. 
since Lx > 0, L ~ 2XLx ~ t^ at most or, (ii) L ~ constant, and 2XLx ~ t^ or faster. 
In the case of possibility (i) the scalar field and in particular the speed of sound plays an 
important role in restricting the small r behavior. In contrast, for possibility (ii) the leading 
behavior is governed by L and the constraints on the absence of super-luminal propagation 
of the scalar fluctuations do not play a role in determining the small r behavior of the metric 
functions. We will discuss each possibility below. 

For possibility (i), it is now straightforward to check that consistency with the Einstein 
equations gives the following behaviors at small r for the metric functions: i^, A ~ r^ or 
faster. At this point the scalar energy density and pressure go like r^ or faster. This behavior 
excludes a purely X independent potential term in the lagrangian for the scalar field since 
that could lead to a constant behavior for small r for p = —L. The sign of the energy 



^To consider an extreme case, if the lagrangian depends only on X then this term is absent. It is then 
trivial to verify that the only solution which is smooth at the origin is ry =constant. This in turn implies 
the constancy of the pressure and the energy density in this case. This is the de Sitter solution. We will not 
discuss this possiblity further. 

^We assume that the lagrangian is Taylor expandable in X. 
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density will now be considered. Equation (TTS!) is useful for this purpose as well. The term in 
brackets on the right hand side of this equation is positive so p is monotonically decreasing 
as we go from some small r to infinity where it approaches the (almost zero) value dictated 
by the near asymptotically flat cosmological solution outside the galactic halo. The pressure 
at small r is therefore positive. Since p is positive and monotonically decreasing to zero for 
all r, one should conclude that because in case (i), p = L — 2XLx ~ t^, the only possibility 
is the trivial solution p = everywhere. While this would be true if the classical theory 
is valid everywhere, it is possible that quantum gravity effects modify the theory at small 
scales negating this argument. The influence of the constraint of subluminal propagation of 
the scalar fluctuations has in fact a very indirect influence on all this. The real reason for 
the absence of the negative energy density is that demanding the leading small r behavior 
of 1]' to be relevant both for the pressure and the lagrangian, eliminates the pure potential 
terms from consideration. We will discuss this in more detail below. 

Possibility (ii) leads to a completely different conclusion. In this case it is easy to check 
that the metric functions have the following behavior for small r: t^, A ~ r^. The essential 
difference with the first possibility is that since eq. flTSl) implies that the maximum positive 
value of p is for small r, therefore, p ~ L > 0. Thus the energy density p = —L can now 
be negative for small r. The lagrangian can contain terms which depend only upon (p and 
not its derivatives. Such potential terms can approach a (negative) constant at small r. 
This case is analogous to the flat space situation discussed in [4j. As we will see below, this 
should not be surprising since the condition for the absence of super-luminal propagation of 
the scalar fluctuations has not played a role here. Even though the energy density can be 
negative in some region, we will now see that under very reasonable assumptions the total 
mass within a large enough region is positive definite. The metric function A is related to 
the mass function m{r) by 

e' = {l- ^-^r. (24) 

r 

Consider the total mass inside a radius R. From ([8]) and (124|) we have, 

i-R j-R 

m = -Att / Lr^dr = -Att (p + 2XLx)r^dr. (25) 

Jo Jo 

In the first term in the integral we perform an integration by parts to write, 

f pr^dr = -- [ "j^r^dr + [-pr^]^. (26) 

Jo o Jq dr 3 

Using (fTSll . and substituting into (!25|) . we obtain for the mass parameter m{R), 

^ = YJ iv'fLxil-^ydr-—\pr']l (27) 

Using (fT3|) the first term is manifestly positive definite and the surface term can be made 
small for large enough R if the pressure falls off faster than -^ (asymptotically near flat 
condition) or if the cosmological solution is such that for large r the pressure is negative as 
is the case with some k-essence models. 

7 



The origin of the negative energy densities in certain regions can be clearly traced to the 
purely potential terms in the lagrangian. Indeed, as goes to a constant for small r, these 
potentials tend to a negative constant. If the total mass parameter is positive, then it 
appears to us that excluding theories where the potential can be negative in some regions is 
not justified. In fact, potentials such as these have been recently considered in a variety of 
situations. For example, in supersymmetric AdS compactifications one encounters potentials 
with local negative maximums, and also a large class of supersymmetric compactifications 
including Calabi-Yau and G2, give rise to effective four dimensional potentials with negative 
regions [H]. The stability of the solutions considered in this paper will be the subject of a 
future investigation. We would like to note, however, that there are many known examples 
where the potential is negative at an extremum and yet the solution is stabilized due to 
gravitational effects [23] as long as the scalar field theory satisfies the positive energy theorem 



Another reason to exclude negative energy densities is its association with super-luminal 
propagation and causality violation and we turn to a discussion of this point. Consider the 
constraints on super-luminal propagation in some more detail. Until now the only constraint 
we have considered is the one arising for the speed of the scalar fluctuations, which is obtained 
from an analysis of the characteristics of the scalar equation of motion. A simple argument 
shows that this has no direct connection with the sign of the energy density. If one changes 
L ^ —L the scalar equations of motion does not change and hence neither does the analysis of 
the characteristics but the energy density changes sign. Therefore there must be an additional 
constraint on super-luminal propagation which is dependent on the sign of the energy density. 
In fact, there have been quite a few papers which investigate the relationship between faster 
than light travel and negative energy densities [TB]. They all differ on the precise definition 
of super-luminal propagation. For the purposes of this paper we follow the discussion of [17], 
which is specific to static spherically symmetric space-times for which the killing time can 
be used to measure the time required for objects or signals to propagate between two of its 
orbits. In [T7| a theorem was proven that if in such a space-time, the (time-like) weak energy 
condition is satisfied, then the signaling time is never faster than the corresponding signal 
in Minkowski space. The normalization of the Killing time is appropriate for an observer 
at very large distances. More specifically, the absence of super-luminal signals in the sense 
defined above, requires T^yt^f^ > for any time-like vector t'^. Demanding the absence of 
super-luminal propagation based on this criterion would eliminate possibility (ii) discussed 
above since it implies negative energy densities at very small r. However, super-luminal 
travel by itself is not threatening as long as there is no causality violation for which closed 
time-like curves are required. In fact, recently [11] it was realized from the viewpoint of 
pure classical field theory, that models which allow for super-luminal propagation even on 
dynamical backgrounds do not necessary posses internal contradictions. In particular these 
models do not lead to any additional causal paradoxes over and above those already present 
in standard general relativity. 

We will therefore take the point of view that neither of the solutions discussed above can 
be a priori discarded (if the total mass parameter is positive definite and well behaved) and 
look for the data to decide if both are realized or not. In this regard it is interesting to note 
that the two possibilities discussed above can predict very different behavior of the rotation 



curves for small r depending on the model for the lagrangian: possibility (ii) implies v^ ~ r, 
and, if such a possibility exists in the complete theory, (i) implies Vc ~ r^ or faster. Our 
discussion suggests that the shape of the rotation curves for small r should therefore provide 
the phenomenological distinction between the possibilities considered. 

If the theories we consider can indeed describe halos of dark matter, then we need to have 
an understanding of its interactions with a black hole which presumably are at the center of 
galaxies. With this in mind, the next section will revisit the "no scalar hair" theorems 
for black holes. 



3 Dressing Black Holes With Scalar Fields 

In this section we study the possibility of a scalar field dressing an asymptotically flat, static, 
spherically symmetric black hole solution in the theory discussed in the earlier sections. As 
discussed in the introduction, this question is of relevance for the stability of the halos 
described previously. 

The distinguishing feature of a black hole is the existence of an event horizon whose position 
at Ts is determined from the condition, g^^{rs) = 0. Physical quantities, in particular the 
components of the energy momentum tensor, are regular at the horizon, from which we 
deduce the regularity of the scalar field and its derivative at r = r^, instead of regularity at 
the origin. With this as the boundary condition, let us consider the scalar field equation of 
motion near the horizon. We want to find the behavior of the allowed solutions consistent 
with the absence of super-luminal propagation of the scalar field fluctuations. From equation 
( TTTj) we see that we need to know the behavior of z/' near the horizon. For this purpose let 
us consider the (tt) component of the Einstein equation ([8]) near r = Vg and write, 

e-^ = AA(r-r,)+0((r-r,)2). (28) 

Thus, 

-A' = h constant + 0{r-rs), (29) 

r — Tg 

and from (^ we obtain that. 



A, = i±i!!i^M > 0. (30) 

rs 

The last inequality follows from the fact that e~'^ must be growing as we move out from the 
horizon. There is the possibility that Ax = 0, which would give the leading behavior in ( l28l) 
as (r — r^)^. This would correspond to the case of an extremal black hole solution, which 
we will not discuss further. We are now able determine the leading behavior of u near the 
horizon consistent with the Einstein equations ([HI [9]). From these, 

(z/' + A') = eVfi;/(r), (31) 

/(r) = -2XLx > 0. (32) 



Integrating this in the neighborhood of the horizon we get, 

v + X = -^f{rs) ln(r - r,) + iT + 0{r - r,), (33) 

where K is an integration constant. Using the previously determined value of A, we finally 
obtain for the leading behavior of v assuming Ax ^ 0, 

u = B ln(r — Tg) + constant + 0{r — r^), (34) 

1 + Krl{L{rs) + fij-s)) 

where B = — -— . (35) 

1 + nrjL[rs) 

It is easy to check from the finiteness of the Ricci scalar at the horizon that 5=1, which 
tells us that /(r^) approaches zero there. Since Lx > 0, we see that X —>■ as r —>■ Vg- The 
exact rate at which this happens will be determined below. From this we see that, 

z/' = h constant + 0{r — Vg), (36) 

r — Tg 

e'' = A,{r - rg) + 0{r - rg)\ (37) 

where the constants in fl36|) and fl29|) are not related and A^ > 0. 

Let us now look for a series solution of the scalar field which is regular at the horizon. In 
addition, we have seen from the discussion immediately preceding Eq. fl36|) . that the solution 
should be such that X = — e^'^0'^ must vanish there. Introducing 6r = r — Vg, we have the 
following expansion: 

ip = 6r' {ao + ai6r + a26r'^ -\ }. (38) 

Substituting this into Eq. fITT]) we get the following indicial equation, 

s{2s + a-3) = 0, (39) 

where, a = ^■^, . . From the discussion on subluminal propagation for the scalar field 
fluctuations in the previous section, a > 1. In this instance, both the solutions for s are not 
immediately excluded provided a is not greater than 3. The special case of a = 3 is the same 
as the s = case except the second solution has a leading log r behavior. Let us first consider 
the possible solutions with near horizon behavior dictated by s = 3/2 — a/2. As a result of 
the smoothness of the action near the horizon, the small br leading order behavior of tl)' is 
br ' ~"' . This behavior is too singular if we are considering only subluminal propagation, 
and hence, this solution is ruled out. We next consider the solution with s = 0. For 
the same reasons discussed in the previous section, ax = 0, and the leading near horizon 
behavior of ip' is of the form 6r. Moreover, this near horizon behavior implies (p' = constant 
or X ~ e~^0'^ ~ 6r. This is completely consistent with the Einstein equation fl?I]) (Recall 
that always, Lx > 0). Indeed from the combined near horizon behavior of u' and A' we 
see that the leading term on the LHS of this equation is a constant, which is the same 
on the RHS. We should comment that the Einstein equation (!3T|) rules out the possibility 
s = 3/2 — a/2 when a < 1 (super-luminal propagation). Thus it should be emphasized that 
consistency (at the horizon) for the s = solution holds irrespective of whether we have 
super-luminal or subluminal propagation of the scalar field fluctuations. 
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We must now check the consistency of the solution at large r. The main question we would 
like to address is whether the solutions at the horizon can match on to the ones at infinity 
in a manner consistent with the Einstein equations. First let us list the behavior of the 
pressure p and -^ at the horizon and asymptotically. We will see how far one can go without 
assuming the weak energy condition. 

Consider, 

JL = -2{- + -)i^Lx, 40 

dr 2 r 

near the horizon. Using ( l36|) and ( l28|) and the fact that (p' = K {a constant), we get, 

^ ~ -AxK^Lx < 0. (41) 

dr 

In order to find its behavior at infinity we need the asymptotic form of z/'. For this purpose, 
one could first find the asymptotic form of A from the integrated version of (jHj): 

e-^ = l-^^^^-^, (42) 

where m{r,rs) = An pr^dr. (43) 

Jts 

Assuming that p falls off faster than -^ at large r, the leading behavior of A' is, 

_^,^ 2Gmir,r^) ^r (44) 

From (13T!) and (132|) it would follow that 

^>'^'"<'-':-' + "'>0. (45) 



which in turn would imply from fj40|) that ^ < also asymptotically for large r. However, 
we want to get information about the large r behavior of ^ without using the asymptotic 
fiatness condition that p falls off faster than -^ at large r. For this purpose we note that 
since ^ < near the horizon, in order for it to be positive asymptotically at large r, it 
would have to vanish in between. We will now show that this is untenable on the basis of 
the Einstein equations, apart from the trivial solution that is a constant. For this purpose 
note that from (l40l) . when ^ = 0, 

- = — . 46 

2 r ^ ^ 

It is easily checked using (Q that this implies, 

1 + Kpr^ 

Since we are looking for regular black hole solutions of the combined scalar-Einstein system, 
we must rule out the possibility of a change in signature of the metric, which means that 
1 + Kpr"^ < or Kpr'^ < — 1. From ( I3T1) and ( !32l) however, 

.2.2 + 3Kpr^ 



r 1 + Kpr"^ 
11 



> 0, (48) 



or, Kpr"^ > —'^. The two constraints on npr"^ are incompatible, hence we conclude that -^ 
cannot vanish between the horizon and infinity and must therefore also be negative asymptot- 
ically. We emphasize that our conclusion that -^ is negative everywhere outside the horizon 
does not rely on either the weak energy condition or the asymptotic fiatness condition. 

Let us next find how the pressure p(r) behaves near the horizon and asymptotically. To 
find p near the horizon, note that here, r/' ~ (r — r^). Thus, at the horizon, p = —p. This 
is the only place we need to worry about whether p is positive or negative. If p is negative 
then p is positive at the horizon and if instead p is positive then p is negative there. Let us 
now consider the possibilities for p at large r. For positive p, since -£; cannot change sign 
in between, p must also be negative at large r. Now it is easy to see that if such a solution 
were to exist, it cannot be asymptotically fiat. Indeed, writing 2?]' Lx = p — L, in (l40l) and 
integrating using the integrating factor we get, 

p = — / {e^r'^Ypdr. (49) 

Let us see what this implies for p{r) at asymptotic values of r. Using the asymptotic condition 
that p falls off faster than ^ at large r, we see that the integral converges and \p\ falls off at 
least like ^ at large r. However, since we have already argued that -£ is negative, it must 
be that p is positive asymptotically at large r. The condition of asymptotic flatness gives 
a positive value for p at large r, thus excluding this possibility. Furthermore, asymptotic 
flatness is incompatible with the constraint Lx > 0. This is a result of the fact that Lx > 
implies that p > —p. Since p is negative in this case and monotonically decreasing, it follows 
from p > —p that p has to be positive and monotonically increasing in order for Lx > to 
be satisfled. The only way we can have an asymptotically flat solution is if p = everywhere. 
We therefore conclude that for positive energy density p, the only allowed black hole solution 
has negative pressure asymptotically at large r, and does not obey the asymptotic flatness 
condition. In k-essence models it is possible to have a cosmological solution in the matter 
dominated epoch with these properties. However, it is not clear how we can match on a 
dynamically evolving scenario onto a static solution. 

When the energy density is negative, the pressure can be both positive or negative for 
large r indicating a solution where one can get either asymptotic flatness or not. However, 
if we invoke the results of [U] regarding super-luminal propagation of signals discussed in 
the previous section then this possibility would be excluded. For reasons mentioned earlier, 
we do not subscribe to this viewpoint as long as the enclosed black hole mass is positive 
and well behaved. We should mention at this point that solutions which can have negative 
energy density somewhere should specially be checked for stability. 

Our discussion in this section suggests that the possibility of a black hole and a scalar halo 
forming at the same time is difficult. This however, does not necessarily make the scalar 
halo model problematic. It is quite possible, and there are other evidence [TS] in support 
of a scenario where a primordial black hole was present before the formation of the galaxy 
itself, triggering the quasar activity. The halo can then be treated as a perturbation on this 
primordial black hole background. An analysis of this is beyond the scope of this paper and 
is the subject of an ongoing project. 
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4 Summary and Conclusions 

We have found some consistent solutions of the scalar-gravity system, which can describe 
galactic halos. Solutions which have negative energy density near the origin have rotation 
curves with f c ~ r at small r. Classically these are the only allowed solutions, however, 
should solutions that have vanishing pressure at the origin be allowed in a complete theory, 
the rotation curves should show much steeper behavior for small r. Solutions with nega- 
tive energy density are associated with super-luminal propagation as discussed in section 
2, however, that does not necessarily imply causality violation. The total energy of these 
configurations can be positive definite even if the energy density is negative somewhere. On 
basis of this, and consistent with similar phenomenon in other physical situations [Hj, we 
have argued against excluding such configurations a priori. Our analysis in section 2 is based 
on spherical halos. However, under reasonable assumptions, we expect similar conclusions in 
general. This is because we can parametrize the departure from sphericity by a parameter 
which should have a smooth limit to zero unless there are topological obstructions. 

We have also reconsidered the question of static spherically symmetric black hole solutions 
in a theory of gravity coupled to scalar fields with non-standard kinetic energy terms. We 
have considered situations where the scalar energy density can be negative in some regions, 
but the total mass is still positive as well as the almost asymptotically fiat and other boundary 
conditions. In the case when the (time-like) weak energy condition is satisfied and the 
asymptotically fiat boundary conditions enforced, we recover the scalar no-hair theorems 
[I9j . Our analysis is based on the minimum of assumptions and different from previous 
ones. We find loopholes not only for the case of negative energy densities, but also when the 
boundary conditions are not asymptotically flat. For reasons mentioned earlier, we do not 
think it is justified to, a priori, discard solutions where the scalar energy density is negative 
somewhere but where the total black hole mass is positive and well behaved. In this context 
it should be mentioned that in [22] consistent black hole solutions have been found in the 
kind of theories considered here which are stationary and do contain regions where the energy 
density can be negative. We have not discussed the stabilty of any possible solutions. This 
and the construction of explicit solutions is also an ongoing project. 



5 Acknowledgements 

We acknowledge Arjun Menon for some initial collaboration and would also like to thank 
him, Fred Adams, Paul Federbush, Saul Teukolsky and Alex Vikman for discussions. We 
are especially grateful to Ted Jacobson for communications about super-luminal propagation 
of the energy flux in general relativity. This work was supported by the US department of 
energy. 



13 



References 

[1] C Armendariz-Picon, V. F. Mukhanov and P. J. Steinhardt, Phys. Rev. Lett. 85, 4438 
(2000) |arXiv:astro-ph/0004134|; C. Armendariz-Picon, V. F. Mukhanov and P. J. Stein- 
hardt, Phys. Rev. D 63, 103510 (2001) | arXiv:astro-ph/0006373j . 

[2] M. Ahshahih a, E. Silverstein and D . Tong, Phys. Rev. D 70, 123505 

(2004) |arXiv:hep-th/0404084|; K Silverstein and D. Tong, Phys. Rev. D 70, 

103505 (2004) |arXiv:hep -th/031022l] ; X. Chen, Phys. Rev. D 71, 063506 (2005) 
|arXiv:hep-th/0 408084j; X . Chen, M. X. Huang, S. Kachru and G. Shiu, JCAP 0701, 
002 (2007) |arX iv:hep-th/ 0605045|. 

[3] A. Sen, Mod. Phys. Lett. A 17, 17 97 (2002) |arXiv:hep-th/0204143 |; A. Sen, 
JHEP 0204, 048 (2002) | arXiv:hep-th/020321l] ; A. Sen, JHEP 0207, 065 
(2002) |arXiv:hep-th/ 0203265|; G. W. Gibbons, Phys. Lett. B 537, 1 (2002) 
farXiv:hep-th/0204008|; G. N. Felder, L. Kofman and A. Starobinsky, JHEP 0209, 
026 (2002) [arXiv:hep-th/0208019j . 

[4] R. Akhoury and C. G. Gauthier. [arXiv:0706.4137[ hep-th] 

[5] N. Ark ani-Hamed, H. C. Cheng, M. A. Luty and S. Mukohyama, JHEP 0405, 074 
(2004). |arXiv:hep-th/0312 099| 

[6] A. Y. Kamenshchik, U. Moschella and V. Pasquier, Phys. Lett. B 511, 265 (2001). 
[arXiv:gr-qc/0103004] 

[7] C. Armendariz-Picon, T. Damour and V. Mukhanov, Phys. Lett. B 458, 209 (1999). 
[arXiv:hep-th/9904075 | 

[8] R. Akhoury and C. G. Gauthier, [iDCivi0804.0420|astro-ph] 

[9] Y. Aharonov, A. Komar and L. Susskind, Phys. Rev. 182, 1400 (1969); R. M. Wald, 
"General Relativity," Chicago, Usa: Univ. Pr. ( 1984)', J- P- Bruneton and G. Esposito- 
Farese, arXiv:0705,4043|gr-qc] . 



[10] C. Armendariz-Picon and E. A. Lim, JCAP 0508, 007 (2005) [arXiv:astro-ph/0 505207| . 

[11] E. Babichev, V. Mukhanov and A. Vikman, JHEP 0802, 101 (2008) 
larXiv:0708.056l( hep-th] . 

[12] For a description of the characteristics of halos see M. Persic, P. Salucci and F. Stel, 
Mon. Not. Roy. Astron. Soc. 281, 27 (1996); P. Salucci, A. Lapi, C. Tonini, G. Gentile, 
I. Yegorova and U. Klein, Mon. Not. Roy. Astron. Soc. 378, 41 (2007) 

[13] T. Matos, F. S. Guzman and D. Nunez, Phys. Rev. D62, 061301 
(2000) [arXiv:astro-ph/0003398|; P. J. Peebles, Astrophys. J. 534, L127 (2000) 
|arXiv:astro-ph/0002495]; L. Arturo Urena-Lopez and A. Liddle, Phys. Rev. D 66, 



083005 (2002) [arXiv: |astro-ph/0207493l ; S. Bharadwaj and S. Kar, Phys. Rev. D 68, 



14 



023516 (2003); J.W. Lee, I. Koh, Phys. Rev. D53 (1996) 2236-2239 ( lhep-ph/9507385| ); 
A. Diez-Tejedor and A. Feinstein, Phys. Rev. D74, 023530 (2006); D. Bertacca, N. Bar- 
tolo and S. Matarrese, JCAP 0805, 005 (2008) a rXiv:0712.0486i 

[14] T. Hertog, G. T. Horowitz, K. Maeda, JHEP, 05, 060 (2003) |hep-th/0304199] . 

[15] K. Scharnhorst, Annalen Phys. 7, 700 (1998) [a rXiv:hep-th/981022l] . 

[16] S. W. Hawking, Phys. Rev. D 46, 603 (1991); K. D. Olum, Phys. Rev. Lett. 81, 

3567 (1998) [arXiv:gr-qc/9805003|; M. Visser, B. Bassett and S. Liberati, Nucl. Phys. 

SuppL 88, 267 (2000) [ arXiv:g r- qc/98 1 0026]; S. Gao and R. M. Wald, Glass. Quant. 
Grav. 17, 4999 (2000) |arXiv:gr-qc7 000702l| . 

[17] B. G. Palmer and D. Marolf. Phys. Rev. D67, 044012 (2003) |arXiv:gr-qc/0211045| . 

[18] R. Bean and J. Magueijo, Phys. Rev. D 66, 063505 (2002) |arXiv:astro-p h/0204486] 

[19] J. D. Bekenstein, Phys. Rev. Lett. 28, 452 (1972); J. D. Bekenstein, Phys. Rev. D 
5, 2403 (1972); M. Huesler, J. Math Phys. 33, 10 (1992); D. Sudarsky, Glass. Quant. 
Grav. 12, 579 (1995); J. D. Bekenstein, Phys. Rev. D 51, 6608 (1995); A. E. Mayo and 
J. D. Bekenstein, Phys. Rev. D 54, 5059 (1996) |arXiv:gr-qc/9602057|. D. Sudarsky 
and J. A. Gonzalez, Phys. Rev. D 67, 024038 (2003) [ arXiv:gr-qc/0207069] . 

[20] K. Lake, Phys. Rev. Lett. 92, 051101 (2004). 

[21] G. W. Misner, K. P. Thorne and J. A. Wheeler, "Gravitation", W. H. Freeman, 
San Francisco, 1973. 

[22] E. Babichev, V. Mukhanov and A. Vikman, JHEP 0609, 061 (2006). 
|arXiv:hep-th/0604075 1 ; E. Babichev, V. Mukhanov and A. Vikman, 
larXiv:0704. 330r[hep-th] . 

[23] W. Boucher, Nucl. Phys. B242 282 (1984) 

[24] Q. H. Park and L. F. Abbott, Phys. Lett. B171 223 (1986) 



15 



